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ABSTRACT 

We compute the circularly polarized signal from atmospheric molecular oxygen. Po- 
larization of Oi rotational lines is caused by Zeeman effect in the Earth magnetic field. 
We evaluate the circularly polarized emission for various sites suitable for CMB mea- 
surements: South Pole and Dome C (Antarctica), Atacama (Chile) and Testa Grigia 
(Italy). An analysis of the polarized signal is presented and discussed in the framework 
of future CMB polarization experiments. We find a typical circularly polarized signal 
(V Stokes parameter) of ~ 50 — 300 fiK at 90 GHz looking at the zenith. Among the 
other sites Atacama shows the lower polarized signal at the zenith. We present maps of 
this signal for the various sites and show typical elevation and azimuth scans. We find 
that Dome C presents the lowest gradient in polarized temperature: ~ 0.3 /iK/° at 
90 GHz. We also study the frequency bands of observation: around v ~ 100 GHz and 
v ~ 160 GHz we find the best conditions because the polarized signal vanishes. Finally 
we evaluate the accuracy of the templates and the signal variability in relation with 
the knowledge and the variability of the Earth magnetic field and the atmospheric 
parameters. 

Key words: atmospheric effects - techniques: polarimetric - cosmic background 
radiation - cosmology: observations 
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1 INTRODUCTION 

The observation of the Cosmic Microwave Background 
(CMB) from the ground can be performed inside the at- 
mospheric windows, far from the main emission lines of 
the most abundant molecules. For the experiments devoted 
to the measurement of the polarization properties of the 
CMB it is mandatory to evaluate, in addition, the effects o f 
the atmosphere. In a recent paper (|Pietranera et all 120071 ) 
the effects of ice crystal clouds in the upper troposphere 
have been considered, but the major effect arises from the 
polarized emission of the Oi molecule. The relevance of 
thi s effect for CMB measureme nts has been pointed out 
by lHananv fc Rosenkranl (|2003h . where for the first time 
a quantitative e valuation of the pola rized signals has been 
presented. Also I Keating et al.l l|l998l ) discuss the relevance 
of the polarized emission of the O2 as a foreground for CMB 
observations. 

Here we compute the polarized emission of the O2 
molecule, develop ing the forma li sm in troduced in the pio- 
neering works by iLenoirl |l967l . [ f968). The interaction of 
the Earth magnetic field with the magnetic dipole moment 
of the atmospheric O2 molecule induces Zeeman splitting 
onto the roto-vibrational lines. As a result of this interac- 



tion, radiation becomes weakly circularly polarized with a 
degree of polarization ~ f0~ 4 — f0~ 5 . A linear polarization 
is also generated, but the signal is three or four orders of 
magnitude lower and still below the current detection sensi- 
tivity. 

This circularly polarized signal can be considered a con- 
taminant for CMB ground-based experiments for two rea- 
sons. The first one is that non-idealities in the instrument 
can produce conversion from circular to linear polarization: 
this leakage from the V Stokes' parameter to {Q,U} can 
be the source of a spurious signal if the circularly polarized 
radiation of atmospheric oxygen is picked-up by the instru- 
ment. Here the effect on the recovery of the cosmological 
polarized signal is evident si nce CMB is expect ed to be only 
(weakly) linearly polarized |Hu fc White!ll997l ). 

The second one is that there is an increasing theo- 
retical effort dedicated to the physical mechanisms able 
to create a non vanishing CMB circular polarization. In 
fact, there are several attempts to add terms to the stan- 
dard electromagnetic lagrangian density, and to predict ob- 
servable effects arising from these extensions. Among oth- 
ers, and without sake of complete ness, we recall: test o f 
the Shiff conjecture on the EEP |Carroll fc Field! 119911 ): 
coupling between photons and slowly varying pseudoscalar 
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fields dHarari fc Sikivid 1199211 | Carlson fc Garretson|[l99l ) 
and (|Finelli fc Galavernil " 2009 ) ; coupling between photons 
and ps eudoscalar fields in p resence of intergalactic magnetic 
fields |Agarwal et alj [20081 ); coupling between photons and 
vector fields to test L orentz-invariance violating processes 
<| Alexander et al.ll2009h . Authors, in all of these cases, show 
how a fundamental process induces optical activity in the 
Intergalactic Medium. In fact, the extensions of the stan- 
dard lagrang ian generate new t erms in Maxwell's equations 
(see eq. 3 in lHarari fc Sikivid (Il992l )). therefore leading to 
modified dispersion relations. We have to mention also the 
possibility of having a V-mode induced by last-sca t tering 
in a magnetized primordial plasma (see iGiovanninil (2009, 
120101 )) and, after recombination, by intergalactic mag netic 
fields through Faraday conversion l|Coorav et al.l l2003). 

It is likely that, with increasing theoretical motivations, 
a dedicated experiment will be proposed in the n ear future to 
updat e the very early upper limits obtained bv lLubin et al.l 
(1983) and iPartridge et all (|l988h . The present work has 
been prepared also in this spirit. 

In section 2 the basic physics of the process is presented, 
while in section 3 the tensor radiative transfer approach used 
to compute the atmospheric integrated signal is described. 
The analysis is applied to different sites: Dome C (Antarc- 
tica), South Pole, Atacama (Chile), Testa Grigia (Italian 
Alps), where mm- wave instruments are located. Results are 
reported in section 4. We compare the effect in different 
places and define the optimal observation bands in order to 
minimize the effect from this contaminant. 



2 POLARIZED SIGNAL THEORY 

O2 is the only abundant atmospheric molecule with non-zero 
magnetic dipole moment, because the two electrons in the 
highest energy state couple with parallel spin. The presence 
of polarized radiation in the atmosphere is associated with 
the interaction of the magnetic dipole moment of the O2 
molecule and the Earth magnetic field. In the millimeter 
wave region of the electromagnetic spectrum, O2 molecule 
presents a set of lines due to roto-vibrational transitions. 
The presence of the Earth magnetic field induces Zeeman 
splitting on these lines, which then become polarized. 

Because of the Zeeman effect a quantum state repre- 
sented by the quantum number J is split in 2J + I states, 
each of them described by a different value of Mj, under the 
condition \Mj\ ^ J. The selection rules permit the follow- 
ing transitions, AJ = +1 (called L + ) and AJ = —1 (called 
L-): 



Table 1 . Line frequency and corresponding quantum number 
used for the computation of the atmospheric O2 polarized signal. 



AJ = +1 
AJ = -1 



J = L 
J = L 



J = L + 1 
J = L-1 



(L + ) 



Under these conditions, selection rules permit to have 
AM j = (tt transitions) and AMj — ±1 (a± transitions). 
In Table [Tj we report the frequencies of O2 roto-vibrational 
transitions allowed by the Fermi Golden Rule and used for 
our calculation, together with their corresponding quantum 
number (I/ ± ). We stop to L = 37 because of lack of in- 
for mation about the li ne broadening and mixing parameters 
(see lLiebe et"aH l|l992l) ) of higher orders. We also verified the 
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effect of adding lines with L > 37: far from the transition 
frequencies the difference is negligible (AV /V < 10 -5 ). 

Radiation intensity and polarization direction are also 
determined by the angle (9) between the line of sight and 
the Earth magnetic field. We must define a reference sys- 
tem onto which the vector E of the electromagnetic wave is 
projected. A good choice is the 3D geomagnetic system (see 
iLenoirl (1 19681 )) centered at the geographic Earth center: z 
axis correspondes to the geomagnetic dipole axis and angles 
ipm and <f> m identify the geomagnetic latitude and longitude 
respectively. Matrix representing the polarized intensity I 
can be evaluated through the coherence matrix (p) for the 
different transitions: 



iAMj (S, L, Mj) = Ptrans,AMj(S, L, M.j) pAMj, 

In this reference frame the matrices for the three different 
transitions become: 



AM j = -1 : p a _ = 
AM j = +1 : p a+ = 
AMj = : pn = 



1 1 cos b 
—i cos 9 cos 2 9 

1 —icos9 
i cos 9 cos 2 9 




sin 2 



(1) 

(2) 
(3) 



where Ptrans denotes the transition probability for the line 
identified by the quantum numbers (S, L, J, Mj) and fix 
their relative intensity. Equations (HJ) - © show that a_ 
and a+ lines produce circular polarization, due to the off- 
diagonal terms in the matrices. Besides the combination of 
the three transitions generate also a small fraction of linear 
polarization. 

The total coherence matrix Atot can be derived by ex- 
pressing the whole contribution given by the superposition 
of the different lines: 
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AMj=+l +J 
btot=C(v,P,T) ]T p AMj jytrans(S,L,Mj,AJ,AMj) 
AMj=—l Mj=-J 

F(v,v k ,Au c ). 

I matrices are multiplied by two weight factors: the first one, 
C(v,P,T), takes into accou nt the depen dence on pressure, 
temperature and frequency (|Lenoir]ll968r ): 



Pv 2 

C(y,T, P) = 0.229— e^; 



(4) 



the second function, F(v, Vk, Av c ), is the line profile. We ne- 
glect the Doppler effect and consider only the collisional line 
broadening, which is dominant in the lower layers of the at- 
mosphere where most of the signal is generated. Therefore 
we u se the collisional line profile (|Van Vleck fc WeisskopJ 
1945t) corrected for t h e line mix ing effects as show n by 
Rosenkranz fc Staelinl (|l988l ) and ILiebe et al l (|l992l ). The 
line profile used is then: 



This matrix expresses the coherence (polarization) proper- 
ties of the signal coming from the Zeeman split O2 lines. We 
get a matrix for a given atmospheric condition of pressure, 
temperature and height. In particular we notice that: 1) 
off-diagonal terms are non vanishing and conjugate purely 
imaginary complex numbers thus radiation is only circularly 
polarized (V 7^ 0) since U parameter is null; 2) circular po- 
larization is larger when the line of sight is aligned with the 
Earth magnetic field: V oc f3 — (a — b) cos 8; 3) the diagonal 
terms are still different (a ~ 7), so that Stokes' parame- 
ter Q oc a — 7 7^ 0; 4) linear polarization is larger when 
the line of sight is orthogonal to the Earth magnetic field: 
Q oc (a + b — c) sin 2 6\ 5) the component of non vanishing 
linear polarization (Q parameter) is nearly aligned with the 
Earth ma gnetic field, even if a small F araday rotation may 
occur (see Rose nkranz fc Stae lin (1988)). Actually it can be 
seen that the ratio o f linear (Q) to circular (V) p olarization 
is Q/V ~ 10" 4 (see lHananv fc Rosenkranz! (|2003l l 1 . In the 
present paper we evaluate only the circularly polarized sig- 
nal, as the linearly polarized one is still below the sensitivity 
of the current experiments (see section 4) . 



F{v,v k ,Av c ) = F A -F Y 



(5) 



where: 

F A = Av c 
F Y = Y k 
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(V k - VY + Al/f (Vk + v) 2 + Avl 
V k —V Vk + V 



K - v y + Avi (v k + v y + Avi 

Here the line broadening is: 



Au c = A P I — 



and the line mixing coefficient is: 



1 , , 300\ „ /300\°' ! 
h = A ■ >•' — J 1 ( — ) 



(6) 



(7) 



Here the coefficients Ap, Ai and A2, fo r each transition, are 
given in table 2 of ILiebe et all l| 19921 ). where calculations 
are compared with measurements. In eq. [5] Uk = vt + Samj , 
where v t is the transition frequency, shown i n table | T] and 
Sam r is a Zeeman induced frequency shift (|Lenoirl 19681 ) 
given a magnetic field B (nT): 



5 am {MHz) = 



f2f (AMj + Mj^ 

K ° B (am j + m.j±±1 



(L + ) 



Here K = 2.8026 x 10~ 5 and AMj = 0, ±1. Taking into 
account all these relations and summarizing them, the total 
coherence matrix can be written as: 



Atot = ap a _ + bp a+ 

a + b i(a — b) cos 9 

-i(a — b) cos 6 (a + b) cos 2 6 + csin 2 9 

a i/3 
-ifi 7 



(8) 



3 COHERENCE BRIGHTNESS 
TEMPERATURE 

3.1 Modeling the atmosphere 

The propagation of radiation emitted by Oxygen molecules 
throughout atmosphere modifies its polarization properties. 
We compute the A matrix via a tensor radiative transfer 
approach. We compute the total coherence matrix in terms 
of brightness temperature and divide the atmosphere in N 
layers of thickness Az — 0.2 km. We interpolate the atmo- 
spheric parameters profiles available in literature to fit this 
thickness. Than we take the average value of the physical 
parameters inside eac h layer a nd us e the theory of radiative 
transfer, described in iLenoirl (|l968T ), to compute the signal 
emerging from one layer, taking into account all the layers 
lying above it: 



,{z,u) =^\NF(Azi,v)Ti. 



(9) 



WF is a weight function for the physical temperature Ti 
of the z-th layer, in which all the information about the 
transfer properties of the layer are enco ded, according to 
the formalism developed in iLenoir! (|l967T ): 

WF(Azi,i/) = P(A%,i/) [l - e - 2A * (,/)A3! *] P*(Azi,u). (10) 

In this relation expresses the total attenuation matrix of 
the i-th layer (see section 2) and P is the product of the 
exponential matrices taking into account the absorption of 
radiation through different layers, that is P = I for i = 1, 
P = e _A ( I/ ) Az » f or i — 2 while for i > 2 we have: 



P(Azi,u) = e 



-A i _ 1 (iv)Az i 



(11) 



To quantify the brightness temperature matrix, we have de- 
veloped a code in which we change the input parameters 
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of the location (geographic co-ordinates) and of the instru- 
ment (pointing directions) considered. The brightness tem- 
perature is therefore computed including the contribution of 
all the transitions of the oxygen molecule, ranging from ~ 
50 GHz to 118.7 GHz shown in Table □ 



3.2 Atmospheric parameters and Magnetic Field 

Vertical profiles for a tmospheric pr e ssure and temperature 



have been taken from 
tion, from |Pe Petrisl 



Tomasi et all (|2006r ) for Dome C star 



_ ( 20081 ) for Testa Grigia, while for At- 
acam43 and South Polqj sites atmospheric parameters are 
available on their related web sites. 

The Earth magnetic field has been calculated with the 
NASA tool based on the IGRF-2010 modefl IGRF is a 
model of the geomagnetic field close to Earth surface, gen- 
erated inside the Earth bodj__. The field can be roughly 
represented as a dipole. Its intensity on the Earth surface is 
~ 3 — 6 x 10 4 nT. IGRF describes the lower spatial frequen- 
cies of the field: IGRF-2010 is truncated at the multipole 
n — 13, corresponding to an angular scale of ~ 15° and to a 
typical wavelength of ~ 3000 km along the surface. Higher 
spatial frequencies, which represent a small scale correction 
to the model, are not accounted for. These components come 
mainly from the magneti zed rocks of the cr ust and are typi- 
cally ~ 5-10 nT rms Csee lHulot et alj (|2002r i). An additional 
contribution comes from the external ring currents present 
in the ionosphere and magnetosphere. The external field in- 
tensity, in quiet conditions of the solar activity, is only a few 
tens of nT. Earth magnetic field is also changing with time. 
The internal field is slowly changing with a global rate of 
~ 80 nT/year. The external field is more rapidly changing 
and increases up to 2 orders of magnitude during magnetic 
storms. Uncertainty of IGRF-2010 model, in a typical quiet 
condition, for the current epoch is estimated to be ~ 10 
nT rms for the total intensity and ~ 1 arcmin for the di- 
rection (see IGRF "health" warning, by F.J. Lowes 4 , IAGA 
Working Group VMOD, January 2010). 

Data concerning geographic, atmospheric and geomag- 
netic parameters of the different observing sites are reported 
in Table [2] Typical atmospheric data and magnetic field are 
evaluated on the ground. In Figure Q] a picture of the lo- 
cal reference frame used to describe the magnetic field B 
direction is shown: p is the angle between B and the local 
horizontal; q is the angle between the horizontal projection 
(H) and the local North. For South Pole X axes is pointing 
towards the geographic longitude £ = 0. 




Y — Eastward 



Z — Downward 



Figure 1. Reference system used to describe the magnetic field 
(B) direction: H is the projection of B onto the horizontal plane. 
The angles q and p define the direction with respect to the geo- 
graphic North (X axis) and the inclination with respect to the lo- 
cal horizontal (orthogonal to Z axis). Downward directions mean 
positive values. 
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Figure 2. Polarized atmospheric O2 signal emitted at the zenith 
for the various observing sites as a function of the frequency. 
The absolute value of V is shown. Colored lines denote negative 
values. Values obtained at Dome C are very similar to those found 
at Testa Grigia, although the sign is the opposite. 



4 RESULTS 

4.1 Frequency templates 

Using the procedure and the code described in the previous 
sections we have been able to evaluate the polarized signal of 
the atmospheric oxygen for each place on the Earth surface, 

1 http: / /www. tuc.nrao.edu/alma/site/Chajnantor/instruments 
/radiosonde/ 

2 ftp : //ftp . cmdl.noaa.gov / ozwv /ozone / Spole / 

3 http:/ /omniweb. gsfc.nasa.gov/vitmo/cgm_ vitmo.html 

4 http://www.ngdc.noaa.gov/IAGA/vmod/igrfhw.html 



assuming that the atmospheric parameters (basically pres- 
sure and temperature profiles) are available. Earth magnetic 
field parameters are computed for the year 2010. 

We computed the circularly polarized contribution (V 
Stokes parameter) for the various sites previously men- 
tioned. The absolute value of the signal at the zenith vs fre- 
quency is shown in Figure [2] The strong contribution of the 
single line at 118 GHz and of the ensemble of lines around 
60 GHz can be immediately noticed. 

Differences among the various sites are due to a number 
of reasons. Among them the most important are: 1. the alti- 
tude of the observing site since the lower layers of the atmo- 
sphere, having an higher pressure (see eq.__]), are the major 
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Table 2. Geographic and atmospheric data on the ground for the various sites considered. For a description of the magnetic field direction 
see Fig. [T] Magnetic field parameters are evaluated for 2010. 



Site 


Latitude (°) 


Longitude (°) 


h (km) 


P (mb) 


T (K) 


B (nT) 


P{°) 


«(°) 


Dome C 


75.06 S 


123.21 E 


3.233 


630 


230 


62713 


-80.8 


-140.4 


South Pole 


90.00 S 





2.800 


680 


242 


55849 


-72.8 


-28.2 


Atacama 


23.01 S 


67.45 W 


5.104 


556 


270 


23296 


-20.2 


-6.3 


Testa Grigia 


45.56 N 


7.48 E 


3.480 


655 


260 


46828 


61.9 


0.7 



contributors to the signal; 2. at mid and low latitudes (Ata- 
cama and Testa Grigia) the amplitude of the Earth magnetic 
field is a factor of two weaker than that at South Pole; 3. 
the amplitude of the effect depends on the scalar product 
between the line of sight and the Earth magnetic field direc- 
tion. For sites close to the magnetic poles the magnetic field 
direction is almost vertical, while at lower latitudes the mag- 
netic field direction is almost horizontal and aligned along 
the North-South direction. 

Black lines in Figure [2] denote values V > (i.e. right 
handed circular polarization), while colored lines denote 
values V < (i.e. left handed circular polarization). We 
can notice that the polarized signal at Testa Grigia has a 
sign reversed respect to other places. This happens because 
the magnetic field B is directed downward in Testa Gri- 
gia (which is in the northern hemisphere), while B is di- 
rected upward elsewhere (all the other sites are located in 
the southern hemisphere). We can notice also a sign rever- 
sal at the frequencies v ~ 99 GHz and v ~ 160 GHz. The 
first null is originated by the 118 GHz line wing crossing 
the 50-70 GHz lines wings. The null at higher frequency is 
originated adding the contribution from lines with L ^ 9. 
Adding lines with higher L the null frequency shifts to lower 
values, converging to the final value shown in Fig. The 
coefficients Aq, A\ and A2 (see eq. 6 and 7), which assume 
different values line by line, seem to play an important role 
in this effect. Anyway, as discussed in section |4~41 it should 
be stressed that above ~120 GHz the model is not yet vali- 
dated, and calculation could be not accurate. 



4.2 Angular templates 

We then built the maps of the polarized signal as seen at 
the various sites in local alt-azimuthal coordinates. Maps of 
the polarized signal, at 90 GHz, seen at the selected sites 
are shown in Figures [3] [4] [5] and [6] Both positive and nega- 
tive values indicate that both right and left handed circular 
polarization is present. We find V = in the maps (red 
contour) at the angular positions where B is orthogonal to 
the line of sight. Elevation cuts along the North-South direc- 
tion and azimuth cuts at constant elevation of 45° of these 
maps are presented respectively in Figures [7] and [8] Observ- 
ing in directions different from the zenith the actual signal 
depends also on the atmospheric thickness. This effect is 
visible in the elevation scans that show a zenith-secant law 
dependence. The spatial structure of these maps is domi- 
nated by the combination of magnetic field direction and 
atmospheric thickness in the line of sight. Each map covers 
a surface, which has a typical scale length of ~ 200 km, at an 
altitude of ~ 30 km, corresponding to a multipole n ~ 200. 
The magnetic field components at these multipoles (or even 
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Figure 3. Map of the atmospheric O2 circularly polarized signal 
at Dome C, at frequency v = 90 GHz. The Lambert projection 
centered at the local zenith, down to 15° of elevation, is shown: 
upward direction is the geographic North. Black circles denote 
elevation intervals of 15°. Color scale ranges from -750 fiK to 
+550 (J.K. White contours denote intervals of 50 (J.K in brightness 
temperature. 



higher) are negligible with respect to the longer wavelength 
components and not larger than the accuracy of the IGRF- 
2010 model (~ 10 nT). 

The effect due to the variation of the atmospheric thick- 
ness is combined in the maps with the direction of the mag- 
netic field. At Dome C, where the magnetic field is almost 
vertical, the two effects compensate each other and the sig- 
nal is flat on a large part of the visible sky. Conversely at 
Atacama, where the magnetic field is almost horizontal, the 
two effects combine increasing the gradient along the visible 
sky. 

The typical signal variation, at 90 GHz, for North-South 
elevation scans (down to 60° from the zenith) is of the order 
of 70 fiK at Dome C, of the order of 300 /J.K at South Pole 
and Testa Grigia, and of the order of 500 [iK at Atacama. 
The typical signal variation for full-360° azimuth scans at 
el = 45° is of the order of 50 fj,K at Dome C, of the order 
of 170 fiK at South Pole and Testa Grigia, and of the order 
of 270 nK at Atacama. 
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Figure 5. Map of the atmospheric O2 circularly polarized signal 
at South Pole. See caption of fig. [3] for details, but notice that 
upward direction is the geographic longitude 1 = 0. 
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Figure 7. Elevation scans of the Oxygen polarized signal for 
the several sites. Cuts are taken along the North-South direction. 
Frequency is v = 90 GHz. 



The signal gradient has also been computed in order to 
give a better description of the signal angular dependence. 
Reference values are summarized in Table [3] The column 3 
of Table [3] presents the minimum and maximum gradients in 
the elevation scans shown in Figure[7] occurring respectively 
at the zenith and at the lowest elevation. We consider only 
the central part where \el\ ^ 45°. It can be noticed that 
the signal variation is minimum (~ 0.3 — 0.8fiK/°) at Dome 
C, while the gradient increases by a factor of 4 at South 
Pole and Testa Grigia and by a factor of 7 at Atacama. The 
column 4 of Table [3] presents the maximum signal gradient 
for the azimuth scans of Figure These values indicate that 
Dome C presents the smoother pattern for the polarized 
emission from atmospheric O2. 
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Table 3. Typical signals and gradients of the Stokes parameter V at the several sites. Signals are evaluated at 90 GHz. Column 2: 
total polarized signal at the zenith (see Fig. [3J|; Column 3: elevation gradient for a North-South scan, minimum (occurring around the 
zenith) and maximum (occurring al lowest elevation) values for el ^ 45° (see Fig. [3; Column 4- maximum azimuth gradient for a scan 
at el = 45° (see Fig.[8]l; Column 5: peak-to-peak polarized signal in a scan at el = 45° (see Fig.[8]l; Column 6: azimuth positions where 
maximum azimuth gradients occur (for azimuth reference see Figs.[3]-[6]|. 



Site 


V x (pK) 




4£ 


SVpp(az)(fj,K) 




Dome C 


-153.4 


0.3 - 0.8 


0.4 


50.2 


129; 310 


South Pole 


-299.0 


1.4 - 3.2 


1.6 


188.1 


61; 242 


Atacama 


-52.5 


2.5 - 5.1 


2.5 


288.0 


83; 264 


Testa Grigia 


+147.1 


1.3 - 2.9 


1.4 


158.8 


90; 271 
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Figure 8. Azimuth scan of the Oxygen polarized signal for the 
several sites. Cuts are taken at constant elevation (45°). Fre- 
quency is v = 90 GHz. 



By means of our template the O2 polarized signal can 
be computed at any time given the site where the observa- 
tions are carried out, giving the possibility of subtracting it 
from the sky signal. For a real experimental setup the sig- 
nals computed for a given direction have to be integrated 
with the efficiency of the bandwidth and convolved with the 
observing beam profile. 



4.3 Bands optimization 

Astrophysical and cosmological signals can be investigated 
by ground based experiments through the atmospheric win- 
dows, where transparency is large. The opacity of the atmo- 
sphere at mm-wave frequency is due to the combined un- 
polarized emission of (mainly) molecular Oxygen (the same 
transitions considered for the polarized emission) and water 
vapor lines (at 22 GHz and 183 GHz). In Fig. [9] the to- 
tal brightness temperature of the atmosphere at Dome C is 
shown between 10 GHz and 200 GHz. 

Inside this interval we define three frequency windows 
of atmospheric transparency: the first between 25 GHz and 
45 GHz, the second between 75 GHz and 110 GHz and the 
third between 125 GHz and 170 GHz. Inside these frequency 
intervals it is possible to observe the sky with a transparency 
larger than 0.8, using a wide bandwidth (Av/v = 0.2 — 
0.3). In these frequency bands we estimate the integrated 
polarized temperature due to atmospheric oxygen. We select 
3 ideal rectangular bandpass which are 10 GHz, 20 GHz and 
30 GHz wide (i.e. corresponding roughly to a 20% bandpass) 




50 100 

Frequency (GHz) 

Figure 9. Brightness temperature of the atmosphere at the 
zenith for Dome C station. The three curves correspond to dif- 
ferent values of the precipitable water vapour (PWV=1 mm cor- 
responding to a typical summer condition, PWV=0.6 mm and 
PWV=0.3 mm corresponding to a typical winter condition). 
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Figure 10. Integrated polarized signal from O2 vs the central 
frequency: bandwidth At/ = 10 GHz; frequency interval from 25 
to 45 GHz; observing towards the zenith. 



respectively for the three bands. The results are shown in 
Figures [TOl 1111 and 1 121 as a function of the central observing 
frequency. 

This analysis allows us to define the optimal frequency 
bandwidths for observation of circular polarization inside 
the windows, while avoiding the largest signals from the 
peaks and the highest part of the line wings. In the 75-110 
GHz interval the high frequency side is preferred because 
here we find a null in the polarized signal. Therefore the 
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Figure 11. Integrated polarized signal from O2 vs the central 
frequency: bandwidth Av = 20 GHz; frequency interval from 75 
to 110 GHz; observing towards the zenith. 
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Figure 12. Integrated polarized signal from O2 vs the central 
frequency: bandwidth Av = 30 GHz; frequency interval from 125 
to 170 GHz; observing towards the zenith. 



optimal frequency band is 87-107 GHz. Inside the 125-170 
GHz interval we do not find a null, but rather a signal that 
becomes smaller and smaller as we rise the central frequency 
of our integration band: for this reason, the contamination is 
smaller in the 140-170 GHz interval. Also in the 25-45 GHz 
interval there is not a null, therefore the farther from the 
line tails the smaller the polarized contribution, and the low 
frequency side is preferred (25-35 GHz). Anyway, it has to 
be kept in mind that a global optimization of the bandwidth 
must also consider the unpolarized signal (see Fig. [9]) and 
the true bandwidth transmission. 



4.4 Accuracy of the templates 

We have estimated the accuracy of the templates at 90 GHz 
for the various sites, taking into account uncertainties and 
typical variability of the principal parameters of the model, 
i.e. the magnetic field and the atmospheric conditions. We 
first considered the magnetic field 4 and summarize the re- 
sults in Table [H taking into account several sources of un- 
certainty or variations: 1) We estimated a typical rms un- 
certainty: <5|B|=20nT, 8p=8q—l'. But at Dome C and South 
Pole we consider 8q=5' . This uncertainty includes both the 
quoted error bars of the IGRF model and the unaccounted 



contribution of the spatial frequencies higher than n = 13. 
Results, obtained using a Monte Carlo for these parame- 
ters, are reported in column 3 of Tabled] 2) Then we evalu- 
ated the secular variation (SV) of the magnetic field, for the 
various sites. We considered the consolidated values for the 
past since 2005 and the forecast values for the next years 
up to 2015. We found a nearly linear variation with time, 
which we reported in column 4 of Table [4] Using the mag- 
netic field values obtained for the several years we computed 
the amplitude of the atmospheric circular polarization. We 
found again a nearly linear variation with time, reported 
in column 5 of Table 2] At Atacama the change of mag- 
netic field direction plays an important role, in particular 
at the zenith, while in all the other sites the strength varia- 
tion is dominating. 3) Finally we evaluated the variation of 
the atmospheric circular polarization in case of an intense 
magnetic storm when the magnetic field strength changes as 
much as ~ 1000 nT. Usually magnetic disturbances last no 
more than a few days and afterwords quiet conditions are es- 
tablished again. Results, which represent unusual maximum 
short-time variations, are reported in column 6 of Table 0] 

We then considered the atmospheric parameters and 
summarize the results in Table [5] We limited our analysis 
to Atacama 1 Dome cQ and South Pole 2 because of the poor 
statistics available for Testa Grigia. We took into account 
several sources of uncertainty or variations: 1) We estimated 
the uncertainty due to the accuracy of the measurements in 
the tempe rature and pressure profile. We took the values 
quoted bv lTomasi et all (|2006l ) (SP exp = 5 mb, ST exp = 0.5 
K). Results (SVexp), obtained using a Monte Carlo distri- 
bution for these parameters, are reported in column 3 of 
Table [5] 2) We then estimated the variability of T and P 
comparing several measured profiles, taken in different days, 
but equivalent conditions. We found a nearly equivalent rms 
variation for the various sites: 8P s t — 0.8 mb, 8T s t — 2.0 K. 
Using these variabilities we obtain 8V 3 t, the signal variation 
reported in column 4 of Table [5] We searched also for a 
day-night variation, but we found that it is limited to the 
temperature of the lowest layers of the atmosphere. The ef- 
fect we found is not larger than SV s t- 3) We also estimated 
the seasonal long-term variability of the atmospheric param- 
eters. We compared measured profiles taken several months 
later and evaluated the maximum variation of T and P. We 
found that at Atacama the variation, even if relevant, is still 
limited to the lowest layers of the atmosphere. Conversely at 
Dome C and South Pole the seasonal variation impacts the 
full air column. In columns 5 and 6 of Table [5] we report the 
typical maximum seasonal variation measured in the lowest 
layers of the atmosphere (the maximum temperature vari- 
ation at the ground is usually even larger). Results on the 
polarized signal (8Vu) are reported in column 7 of Table [S] 

As shown by Tables [3] and [S]the accuracy is of the order 
of fiK. Larger effects arise in case of large magnetic distur- 
bances or seasonal changes. We did not take into account 
the line profile model. The line frequencies are known with 
negligible error bars, but the Oxigen absorption values have 
error bars of at least 5% at 90 GHz. This accuracy applies up 
to ~ 120 GHz, but at higher frequencies model and values 
of Oxigen absorption are not yet validated and calculations 
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Table 4. Typical polarized signal uncertainties or variations due to the amplitude od the Earth magnetic field. We evaluated the effect 
at the zenith and at (e(=45° ; az=0). Column 3: rms error bar due to the IGRF model uncertainties. Column 4- Secular variation of 
the magnetic field. Column 5: change of the signal due to SV of magnetic field. Column 6: signal variation due to a typical big magnetic 
disturbance (md). 



Site 


Elevation 


SVrms(fJ,K) 






SV md (jj,K) 


Dome C 


90° 


0.05 


-21 


+0.05 


~ 3 


Dome C 


45° 


0.07 


-21 


+0.13 


~ 2 


South Pole 


90° 


0.11 


-63 


+0.38 


~ 5 


South Pole 


45° 


0.16 


-63 


+0.44 


~ 7 


Atacama 


90° 


0.06 


-52 


-0.27 


~ 2 


Atacama 


45° 


0.17 


-52 


+0.23 


~ 9 


Testa Grigia 


90° 


0.07 


+26 


+0.12 


~ 3 


Testa Grigia 


45° 


0.07 


+26 


+0.04 


~ 1 



Table 5. Typical polarized signal uncertainties or variations due to the atmospheric parameters. We evaluated the effect at the zenith and 
at (ei=45° ; az=0). Column 3: rms error bar due to temperature and pressure profile measurement uncertainties (<5P = 5 mb, ST = 0.5 
K). Column 4- signal variation due to short-term rms variation of temperature and pressure. Column 5: typical seasonal (long-term) 
temperature variation. Column 6: typical seasonal (long-term) pressure variation. Column 7: signal variation due to a typical seasonal 
(long-term) variation. 



Site 


Elevation 


SV exp (fj,K) 


SVst^K) 


ST lt (K) 


SP u {mh) 


SV lt (jtK) 


Dome C 


90° 


0.7 


0.5 


36 


15 


3.4 


Dome C 


45° 


0.7 


0.5 


36 


15 


2.9 


South Pole 


90° 


1.1 


0.8 


26 


15 


10.0 


South Pole 


45° 


1.4 


1.1 


26 


15 


13.7 


Atacama 


90° 


0.1 


0.1 


10.5 


2.6 


0.23 


Atacama 


45° 


0.6 


0.3 


10.5 


2.6 


0.9 



should be used with caution. The effect of these uncertain- 
ties should be added to the accuracy quoted previously. It 
is also important to stress that the computed uncertainty 
refers to the total power polarized signal. Usually we are 
interested to subtract signal gradients between different ob- 
serving directions, and the related uncertainty is definitely 
smaller than the total estimated value of SV, in particular 
for small angular scales. 



5 CONCLUSIONS 

In this paper we evaluate the contribution of atmospheric 
O2 Zeeman split roto- vibrational lines in the mm-wave part 
of the electromagnetic spectrum, seen as a contaminant for 
CMB polarization experiments. A code to simulate this sig- 
nal for various sites (Dome C and South Pole in the Antarc- 
tica, Testa Grigia in the Italian Alps and Atacama, Chile) of 
interest for mm-wave astronomy has been deve loped on the 
basis of the theory developed in iLenoirl l|l967t l and iLenoirl 
(1968)). In this way we compute a template of the atmo- 
spheric circularly polarized signal for every site provided 
that the vertical profiles of atmospheric pressure and tem- 
perature are available. 

Maps of the circularly polarized signal in local coor- 
dinates have been obtained for each site examined. These 
maps show a dipole-like large scale feature with a maximum 
gradient aligned approximately, for low-mid latitude sites, 
with the North-South direction, and for polar sites, along the 
magnetic field lines. Atacama is the examined place where 
the polarized signal at the zenith is the lowest. This is due to 



the altitude of the site and to the amplitude and direction 
of the magnetic field. Unfortunately scanning experiments 
at Atacama suffer a larger polarized gradient with respect 
to other sites. In fact the angular variation of the signal is of 
order ~ 1.3 — 5.1 ^K/° at 90 GHz for all sites except Dome 
C, where it is a factor ~ 4—7 lower, because here the atmo- 
spheric polarized signal is the flattest, along the full local 
visible sky. 

We then evaluated the integrated polarized signal inside 
an ideal bandpass of a wide-band (Au/v = 0.2) experiment 
for 3 mm-wave atmospheric windows between 10 and 200 
GHz. We found the frequencies where it is more convenient 
to observe, because here the signal vanishes or decreases 
significantly. Finally, an estimate of the template accuracy 
and variation has been done, using the main features of the 
magnetic field and the statistics of measured atmospheric 
profiles. At 90 GHz, we obtained a typical uncertainty of 
the order of ~ 1 fiK, dominated by accuracy and variability 
of the atmospheric parameters. Large variations (up to ~ 10- 
15 fiK) are expected in case of large magnetic disturbances 
and due to seasonal atmospheric variations. Uncertainty and 
secular variation of the magnetic field have, on the other 
hand, a negligible effect. 

A linear polarization is also generated and it is larger 
when the line of sight is orthogonal to the Earth magnetic 
field. However the ratio of linear to circular polarization is ~ 
10 -4 , therefore still below the current detection sensitivity. 
Moreover the non vanishing component of linear polarization 
is nearly aligned with the Earth magnetic field. 
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